Linear superpositions of gap solitons in periodic Kerr media 
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The existence of a novel type of solitons in periodic Kerr media constructed as superposition of nonintcracting 
gap-solitons of different kinds (bright, dark and periodic) is first demonstrated. The periodic modulation of the 
nonlinearity is used to suppress the cross phase modulation between components to allow the superimposed 
beam to propagate for long distances. © 2011 Optical Society of America 
OCIS codes: 190.0190, 190.6135. 



Photonic structures with periodic modulations of pa- 
rameters are ideal systems to investigate nonlinear prop- 
erties of wave dynamics (see e.g. [I]). In this context, ar- 
rays of optical waveguides are one of the most used con- 
figurations for practical implementations of various phys- 
ical phenomena involving optical solitons. To avoid prob- 
lems of cross-phase modulations affecting beams con- 
sisting of different modes, which differ by frequencies 
and/or by wavelengths, the nonlinear properties of the 
light propagation in these structures have been inves- 
tigated mainly for monochromatic light. Instability in- 
duced by cross-phase modulation, however, could arise 
also for monochromatic light, if linear superposition of 
nonlinear modes of different wavelengths are considered. 
Meantime, stabilizing effects on the dynamics achieved 
by means of periodic modulations of the nonlinearity 
have been explored in the context of Bosc-Einstein con- 
densates in optical lattices where it has been shown that 
periodic modulations of the interatomic interaction allow 
one to induce Bloch oscillations [2], intra-band (Rabi) 
oscillations j3j and dynamical localization [4] of gap soli- 
tons which persist for long times (long propagation dis- 
tances in the optical context). In view of the analogies 
between photonic lattices and Bose-Einstein condensates 
in optical lattices it is of interest to investigate the pos- 
sibility of using periodic modulations of the nonlinearity 
as a tool for controlling the cross phase modulation in 
superimposed beams. 

The aim of this Letter is to demonstrate the existence 
of a novel type of solitons in periodic Kerr media made 
as superpositions of noninteracting gap-solitons of dif- 
ferent wavevectors and different kinds (bright, dark and 
periodic). The periodic modulation of the nonlinearity 
is used to suppress the cross phase modulation between 
components to allow the superimposed beam to propa- 
gate for long distances. 

As model equation for the propagation of a monochro- 
matic beam along z-axis in a periodic Kerr media we use 
the following nonlinear Schrodingcr equation for the di- 



mensionless electric field envelope q (see e.g. [T]): 

iQz = -qxx + V(x)q + G(x)\q\ 2 q. (1) 

Here V(x) = V(x + it) and G(x) = G(x + it) denoting 
the linear refractive index and the nonlinearity modula- 
tions in the transverse direction, respectively. Hereafter 
we fix V(x) — T^cos(2x) and G{x) = a + Gcos(2x), 
both taken of period it without loss of generality. The 
average Kerr index a, as well as amplitudes of lin- 
ear refractive index, V, and of nonlinearity, G, modu- 
lations are considered as free parameters, used below 
for achieving zero cross phase modulation conditions. 
The periodicity of the linear refractive index induces a 
band structure determined by the eigenvalue problem 
d 2 tp a ,k/dx 2 — V(x)<p a ,k = ba(k)<p a ,k, where the eigen- 
values b a (k) (propagation constants) are periodic func- 
tions of the Bloch vector k G [—1,1] and the associated 
orthonormal set of eigenmodes (p a ,k satisfy the Bloch 
condition (p a .k( x + ?r) — e <p a ,k(%) (here a denotes 
the band index). We are particularly interested in the 
modes bordering different edges of the same or of differ- 
ent bands. Respectively, we use subscripts /, u to denote 
lower and upper propagation constant values (bi < b u ). 
The self-phase modulation of small amplitude soli- 
tons bifurcating from band edges &z jU is described by 
Xi.u = Jo G(x)\ipi tU \ 4: dx, while the cross-phase modu- 
lation is determined by \ = f£ G(x)\<pi\ 2 \<p u \ 2 dx. Then 
necessary conditions for the existence of the bright soli- 
tons, combined with the requirement for zero cross-phase 
modulation, allowing them to bifurcate independently 
from different edges, take the form 



Dixi < 0, D uXu < 0, % = 0, 



(2) 



where Di M denote the linear diffraction coefficients of 
the respective I, u modes: Di_ u = —d 2 b^ u /dk 2 (see [2] 
for details). These conditions can be satisfied for proper 
choices of V(x) and G(x). 

Let us now consider a superposition of two localized 
solutions, wi^{x)e Lhl - 2Z of Eq. (fTJ), i.e. 

q(x, z) = w 1 (x)e lblZ + w 2 (x)e lb2Z+ie . 



(3) 
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Fig. 1. (color on line) (a) Families of the modes bifurcat- 
ing from the lower, U\, and upper U 2 edges of the first 
gap: b 6 [—2.1659,0.7332] (the edges coincide with the 
panel boundaries) ; (b) U\ vs U 2 for the solitons fulfilling 
(|4]); (c— f) Intensity distributions \q{x, z)\ 2 and (g-j) spec- 
tra \Q(b, k)\ 2 at distance zq = 1000 of input beams © 
consisting of two superimposed solitons with parameters 
(c,g) Ux = 0.32 (&! = -2.15), U 2 = 0.143 (b 2 = 0.71) 
(d,h) C/i = 0.645 (fo = -2.1), C/ 2 = 0.33 (b 2 = 0.6) 
(c,i) t/i = 2.77 (6i = 0.732), C/ 2 = 1.7 (b 2 = 0.938) 
(f,j) [/i = 2.77 (&i = 0.732), U 2 = 5.28 (& 2 = 1.0). In 
panels (c,d,g,h) solitons are near edges of the first gap 
while in panels (e,f,i,j) near edges of the first band. Other 
parameters are a = — 1, V = —3, G = 3.05 (a-d,g,h), 
G = 1.282 (e,f,i,j). 



Here without loss of generality the functions Wj (x) (j = 
1, 2) are taken to be real solutions of the stationary equa- 
tion d 2 Wj/dx 2 + bjWj — V(x)wj + G(x)Wj = while 
the constant O denotes a relative phase. For propaga- 
tion constants sufficiently close to the gap edges, i.e. 
for b\ w bi and b\ « b u , these solutions can be rep- 
resented as wi i2 (x) = Ai, 2 (ni, 2 x)tpi tU (x), where /xi, 2 = 
\/\bi. 2 — bi, u \ <C 1 is a small parameter characterizing 
detunings to the gap, and A,(-) arc slowly varying am- 
plitudes. Then ^j ensures that \ ~ A*j *C Xl,u so that 
the cross phase modulation of Bloch states implies strong 
suppression of the cross-phase modulation of the small 
amplitude solitons Wi j2 (x). Respectively, we expect the 
superposition of such solutions to approach a stable so- 
lutions of the nonlinear equation. When the detuning 
towards the gap increases, the condition @ must be 



modified, to ensure 

/oo 
G(x)w 2 w 2 dx = 0. (4) 
-oo 

To show that when (0} are satisfied the su- 
perposition ([3]) becomes highly accurate approxima- 
tion of the exact solutions we recall that Eq. (TTJ) 
admits two conserved quantities: the energy flow 

U[q] 



J_ , Iqpeta; and the Hamiltonian H[q] 



J-oo [\<lx\ 2 + V(x)\q\ 2 + \G{x)\q\ A ] dx. For the solutions 
j4{, then the above conditions imply: U[q] — U[w\] + 
U[w 2 ] and H[q] = H[u>i] + H[w 2 ] so that, in spite of the 
nonlinearity in Eq. (TTJ), the superimposed beam ([3]) will 
propagate as an exact solution, without appreciable dis- 
tortions. In the case of solitons with 6i i2 bifurcating from 
opposite edges of the first gap [Fig.QJa)], the first of con- 
ditions (@| is satisfied due to the parity of solutions W\ (x) 
(odd) and w 2 (x) (even). As to the cross-phase modula- 
tion, for fixed values of V and G, for each even soliton 
with b 2 (and the energy flux U 2 — U[w 2 ]) there exists 
its unique odd counterpart with b\ (and the energy flux 
U\ = U[wi\), satisfying the second condition in ^. We 
remark that for two low-amplitude solitons the second 
condition in (|4]) is usually satisfied for G > Go , where Go 
denotes the amplitude of G(x) for which \ = (for the 
refractive index parameters of FigJTJa) Go ~ 2.98 while 
G = 3.05). For this situation the respective dependence 
U\(U 2 ) is illustrated in Fig.[IJb). We however emphasize 
that only stable points of the presented curve can be 
used in an experiment. Typically, only small- amplitude 
solitons are stable: U\ < 0.5 for the parameters chosen 
in Fig. ffl 

The fact that the superposition can propagates with- 
out significant distortion is confirmed by Figs. [Uc,e). 
In order to emphasize the distortionless propagation of 
the beam, at z = 1000 (denoted in Figs. QJc-f) by 
horizontal lines) we applied a refractive index gradi- 
ent in the transverse direction, modeled by adding the 
term O.OOlxq^a;, z) in the right-hand side of Eq. (TTJ) 
for z > zq. We observe that when conditions ((4]) arc 
met [Figs. [ljc,e)] the two superimposed solitons propa- 
gate undistorted until the refractive gradient is applied 
and after that it splits into two gap solitons propagat- 
ing in opposite directions, due to the opposite signs of 
the diffraction coefficients Di, u . If, however, Eq. (|4|) is 
strongly violated we get that the superimposed beam 
is destroyed after relatively short propagation distances 
(i.e. for z ^C Zq) [see Figs. [TJd,f)]. The spectra Q(b, k) = 
/_ l( x t z)e~ tbz ip a ^(x)dzdx of the superimposed beams 
computed at z — zq are shown in Figs.[TJg-j) (each panel 
corresponding to the spatial dynamics shown in the mid- 
dle row just above it). We see that when the zero cross 
phase modulation condition is satisfied [Figs. [T{g,i)] the 
output spectrum is centered around the same parameters 
b 1 k in the Fourier space as the respective input beam 
(in spite of the very long propagation distance), while in 
the opposite case, [Figs. Qlhj)], the spectrum becomes 
strongly deformed. Note the differences between spectral 
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Fig. 2. (color on line) Propagation of two superimposed 
gap solitons with centers shifted at the input by Ax = 
4-7T (panels a,c), or by Aa; = 107r (panels b,d). Soliton 
parameters are the same as in Fig. [He) (panels a,b) , or 
FigOJe) (panels c,d). 



widths and transverse spatial extension of the superim- 
posed beam in Figs. [ljc,g) and the ones in Figs. [TJc,i)- 
While for the former case (component solitons at oppo- 
site edges of the first gap) the orthogonality condition is 
always satisfied due to the parity, in the latter case (com- 
ponent solitons at opposite edges of the first band) this 
is true only in the small amplitude limit, as we remarked 
before (in this case the orthogonality of the solitons can 
be achieved only through the orthogonality of the Bloch 
states). This explains the wider width in the transverse 
direction of the latter beam and the narrowing of the 
spectral widths of the component gaps-solitons. 

In Fig. [2] we have shown the effect of an initial dis- 
placement Ax on the propagation of the two superim- 
posed solitons in Fig. QJc,e). We see that the super- 
positions remain stable both for small [Figs|2la,c)] and 
large [Figs|2fb,d)] displacements. This is understood by 
the fact that, being the solitons of small amplitudes, 
nonlinear-orthogonality condition Q is only slightly vi- 
olated by small displacements. 

The above results can be extended to superpositions of 
nonlinear periodic waves with a bright soliton [Fig. |3Ja- 
c)] and to superpositions of dark solitons [Fig.[3Jd f)]. In 
the first case the periodic (even) wave w^{x) bifurcates 
from the lower edge of semi- infinite gap [Fig.|3][a)], while 
the bright (odd) soliton w\ (x) bifurcates from the lower 
edge of first gap [Fig. [3th)]. The both entities are stable; 
we however note the difference in their amplitudes [c.f. 
the scales of \q\ 2 in Figs. [3][a) and[3fb)]. In the second 
example, the dark solitons bifurcate from the opposite 
edges of the first gap towards the first [Fig.JHJd)] and the 
second [Fig.[3][e)] allowed bands. In both cases the condi- 
tion for zero cross phase modulation is the same as in Q 
(or x — in the small amplitude limit) but conditions 
for existence and stability of the modes must be changed 
as Dixi < 0, D u Xu > and A,«Xi,« > ° for periodic 
wave-bright soliton and dark soliton - dark soliton su- 
perpositions, respectively. In the bottom panels of Fig. [3] 
are shown the output spectra computed at z = 1000 of 
the corresponding superimposed modes. In both cases 
we see that the linear superposition propagates in sta- 
ble manner for a long distance as consequence of the 
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Fig. 3. (color on line) (a,b) Even periodic background 
with the energy flux per period C/ 2 = $ Q w\{x)dx — 
0.103 and b 2 — 0.76 (panel a), superimposed with an 
odd soliton having U\ = 0.132 and b\ = —2.16 (panel 
b) for the parameters a = — 1, V = —3, G = 6; (d,e) 
Superposition of two dark solitons bordering edges of 
the gap b G [—2.1659,0.7332]: even dark soliton with 
b 2 = 0.78 (panel d) and odd dark soliton with b\ = —2.18 
(panel e) for V = 3, a = 1, G = 2.95; (c,f) Spectrum 
\Q(b, k)\ 2 relative to superpositions of the waves in (a,b) 
and in (d,e) respectively, both taken at z = 1000. 



zero cross phase modulation between beam components. 
Similar results can be derived for solitons subjected to 
refractive indexes and nonlinearity modulations of dif- 
ferent periods, as the ones considered in Ref. [5]. 

In conclusion, we have demonstrated the possibility of 
stable propagation of new types of superimposed beams 
in periodic Kerr media which are constructed as linear 
superpositions of gap solitons of several types and dif- 
ferent wavelengths. The Kerr nonlinearity was designed 
in such a manner that the cross-phase modulation be- 
tween beam components is eliminated or strongly re- 
duced. These results can be of practical interest for 
multi-component parallel transmission and manipulation 
of optical signals in nonlinear media. 
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